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Abstract — We demonstrate numerically the diffraction- free propagation of sub-wavelength 
sized optical beams through simple elements built of metal-dielectric multilayers. The proposed 
metamaterial consists of silver and a high refractive index dielectric, and is designed using the 
effective medium theory as strongly anisotropic and impedance matched to air. Further it is 
P-H . characterised with the transfer matrix method, and investigated with FDTD. The diffraction- 

's^ ■ free behaviour is verified by the analysis of FWHM of PSF in the function of the number of 

periods. Small reflections, small attenuation, and reduced Fabry Perot resonances make it a 
flexible diffraction-free material for arbitrarily shaped optical planar elements with sizes of the 
q 1 order of one wavelength. 

Oh 

1. INTRODUCTION 

Since the seminal paper by Pendry pQ, subwavelength imaging at visible wavelengths has been 
demonstrated in much thicker low-loss layered silver-dielectric periodic structures [21 13 SI El El El 
ElE]. In fact, the effective medium theory (EMT) provides sufficient means to explain enhanced 
transmission through the metal-dielectric stacks. Subwavelength resolution results from the extreme 
anisotropy of the effective permittivity tensor [21E1E]- 

Let us continue the introduction by linking the metal-dielectric multilayers for sub-wavelength 
imaging with the concepts taken from Fourier Optics. We refer to the model of a linear shift 
invariant scalar system (LSI) [TOl [TT] for the description of optical multilayers in a situation when 
they act as imaging nano-elements for coherent monochromatic light. In-plane imaging through a 
\& ' layered structure consisting of uniform and isotropic materials represents a LSI, for either TE or 

TM polarisations. Linearity of the system is the consequence of linearity of materials and validity 
(S| ■ of the superposition principle for optical fields. Shift invariance results from the assumed infinite 

\ perpendicular size of the multilayer and the freedom in the choice of an optical axis. 

A scalar description is valid for the TE and TM polarisations in 2D since all other field compo- 
nents may be derived from E y or H y , respectively, where the co-ordinate system is oriented as in 

>: fig.m 

For the TM polarization, the magnetic field H y (x,z) may be represented with its spatial spec- 
trum H y (k x , z) 

H y (x,z)= / H y (k x ,z)exp(ik x x)dk x , (1) 

J — oo 

where, at least for lossless materials, the spatial spectrum is clearly separated into the propagating 
part k% < k^e and evanescent part k^ > k^e. The transfer function t(k x ) (TF) is the ratio of the 
output to incident fields spatial spectra and corresponds to the amplitude transmission coefficient 
of the multilayer 

H y (k x , z = L) = t(k x ) • Hy nc {k x , z = 0). (2) 

Due to reflections, the incident field Hy nc (k x , 2 = 0) differs from the total field H y (k x , z = 0). 

The point spread function (PSF) is the inverse Fourier transform of the TF and has the inter- 
pretation of the response of the system to a point signal S(x). The response to an arbitrary input 
Hy nc (x, z = 0) can be further expressed as its convolution with the PSF 

H y (x, z = L) = H 1 ™^, z = 0) * PSF(x) (3) 
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Figure 1: Schematic of a periodic silver-dielectric multilayer with symmetric (a) or non-symmetric (b) 
composition, and the equivalent effective medium slab (c). 
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Figure 2: a Intensity transmission and reflection coefficients (T, R) of the symmetric multilayer as a function 
of the filling factor, plotted over the transfer function (in vertical cross-sections of the color map). Phase 
isolines are separated by 7r/4. The pitch A = d\ + cfe is fixed at 22nm, and the total thickness is L = 30 A. 
b. Characteristics of the symmetric and non-symmetric multilayers, and of the equivalent slab made of an 
effective medium, in the function of the number of periods. Intensity transmission and reflection coefficients 
calculated at normal incidence (top); FWHM of PSF and \PSF\ 2 (bottom). 



PSF of an imaging system usually provides clear information about the resolution, loss or en- 
hancement of contrast, as well as the characteristics of image distortions. However, for subwave- 
length imaging, the PSF is not a straightforward measure of resolution and even imaging of objects 
smaller than the FWHM of PSF is possible P2j. 

In this paper we demonstrate a diffraction-free material for subwavelength sized optical beams. 
We combine the following properties of the multilayer: PSF with sub-wavelength size and little de- 
pendence on the thickness of the structure, high transmission, low losses, and a limited dependence 
of the imaging properties on the size of external layers. Together, these properties allow to use a 
multilayer as a flexible construction material for various optical imaging nano-devices. 

2. IMAGING WITH SUB- WAVELENGTH RESOLUTION IN METAL-DIELECTRIC 
MULTILAYERS 

The dispersion relation of a two-component infinite stack for the TM polarisation has the form [13] 
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( if a "\ (i *k\ (i n -p\ w sm(k lz fA)sm(k 2z (l-f)A)fk lz e 2 k 2z e 1 \ 
cos(i^A) = cos(k lz fA) cos(k 2z (l - f)A) 1^- + —j^- J , (4) 

where K z is the Bloch wavenumber, A = d\ + d 2 is the period of the stack, d{ and q are the layer 
thickness and permittivity of material z = 1,2, the filing fraction of material 1 is / = d\j A and the 
local dispersion relations are kf z + k\ — k^ei. Wavevector k x is conserved at the layer boundaries 
and depends on the incidence conditions. The Bloch wavenumber K z is real for a Bloch mode in 
a lossless stack, however for evanescent waves in a finite stack or for lossy materials K z may be 
complex. In the first BZ, the real part of K z satisfies tt/A <= re(K z ) <= tt/A. The group velocity 
as well as the imaginary part responsible for absorption do not depend on the choice of BZ. 

When the layers are thin K Z A, k{ z A « 1, the second order expansion of (jlj) over the arguments 
of trigonometric functions leads to the dispersion relation for an uniaxially anisotropic effective 
medium 

K 2 Je x + k 2 Je z = kl (5) 

with e x = fei + (1 — f)e 2 and e z = l/(f /e\ + (1 — f)/e 2 ). This is the basis of the effective medium 
approximation thoroughly discussed by Wood et al. [2], also as the basis for applying the near- 
field approximation. The transmission coefficient of the Fabry-Perot (FP) slab consisting of the 
homogenized effective medium is then given as [2], 

t(k x ) = (cos(k z L) - 0.5i(K z e x /k z e + k z e/K z e x ) sin(^L)))- 1 , (6) 

where k z and e refer to the external medium, and L is the total thickness of the slab. At the 
same time t(k x ) is the already mentioned coherent amplitude transfer function of the imaging 
system [IT]. When \e x /e z \ « 1 (which is satisfied for a lossless metal for ei/e 2 = —d\jd 2 [7J), for 
certain slab thickness the resonant FP condition becomes independent on the angle of incidence. 
This happens when L = Am/^m E N. Then t{k x ) = exp{—iK z L) and the FP slab introduces 
the same phase shift for all harmonics of the spatial spectrum. Belov and Hao [7J proposed to 
combine this condition with impedance matching between the external medium and the effective 
FP slab e — e x — fe\ + (1 — f)e 2 and referred to that regime as canalization. However, Li et al. [5] 
questioned the importance of impedance matching in favour of the FP resonance condition. In fact, 
for a lossless metal and dielectric, the FP resonance is sufficient to entirely eliminate reflections 
resulting in perfect imaging t(k x ) = 1 without impedance matching. In reality, losses make the 
condition e — e x only approximate and, at the same time, the finite value of A limits the validity 
of homogenisation. Moreover, transmission through a finite slab strongly depends on the material 
and thickness of the external layers and appears to be the largest for a symmetrically designed 
slab (Fig. [Jji) with half- width dielectric layers located at the boundaries [3j. 

3. NUMERICAL DEMONSTRATION OF DIFFRACTION-FREE PROPAGATION OF 
SUB- WAVELENGTH SIZED OPTICAL BEAMS 

After recalling the theory of sub- wavelength imaging in metal-dielectric multilayers, we now focus 
on an imaging regime which may be called diffract ion- free. We note, that the FP resonances are 
accompanied with a field pattern inside the slab similar to a standing wave [21 [9] . Therefore look 
for structures for which t(k x ) ~ const nonetheless FP resonances are weak. 

Let us now focus on an example of a metal-dielectric periodic multilayer consisting of silver and 
high refractive index dielectric, which enables us to demonstrate and explain the diffraction-free 
propagation of sub-wavelength sized optical beams. The structure operates at the wavelength of 
A = 422nm, when the permittivity of silver is equal to t\ — — 5.637+0.214z[14], and the permittivity 
of the dielectric such as Ti0 2 or SrTiO% jTUl ES] is e 2 = (2.6) 2 . Layer thicknesses are assumed 
to be d\ — lOnm, d 2 — Ylnm with periodic symmetric or non-symmetric composition shown in 
fig. []Ji,b. The fraction d\jd 2 may be justified with the use of EMT. The corresponding permittivity 
of the effective medium gives c x 1.02 + 0.1 z and c z ~ — 158 + 191z, which assure impedance 
matching with air e = 1 « e x together with the condition for the extreme anisotropy \e x /e z \ << 1. 
The imaginary part of the refractive index n x — is reduced by the factor of 50 compared to 
bulk silver, resulting in low-loss transmission. In flgE^t we show the transfer function calculated 
rigorously with TMM for a range of filling factors < d\/A < 0.6, when the total thickness 
of the structure is fixed at L = 660nm. Indeed, the impedance matching which occures when 
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Figure 3: Diffraction- free transmission of the subwavelength sized beam through the multilayer (FDTD 
simulations, with the incident CW beam limited by a subwavelength aperture in a perfect conductor): 
a. time-averaged Poynting vector S z in a rectangular slab; b,c. instantaneous magnetic field H y in a slab 
with layer boundaries oriented at 30° towards the external boundaries for various aperture sizes. 



Instantanous distribution of Hy 



Instantanous distribution of Hy 



-1.5 



-2.5 




n 



0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 



0.5 1 1.5 2 2.5 3 Hy 

z IX [a.u.] 



n 



a. 



b. 



0.5 1 1.5 2 2.5 Hy 

z IX [a.u.] 



Figure 4: Simple optical nano-elements internally made of the diffraction- free material a. double slab; b. 
prism;. 



d\ — lOnm, d,2 = Vlnvn results in reduced reflections, high transmission of both propagating and 
evanescent spatial harmonics, and a flat phase of the transfer function for a broad range of k x /ko. 
The size of the corresponding PSF is of the order of A/ 10 « 40nm assuring the super-resolving 
properties of the device. Fig. Eh shows the FWHM of PSF alongside with the reflection and 
transmission coefficient in the function of number of periods, for the symmetric, non-symmetric and 
effective- index composition of the structure. The major properties of the structure are the following: 
the size of PSF varies slowly with the size of the structure and is almost the same for the symmetric 
and non-symmetric composition, the FP resonances observed in transmission are weak (as opposed 
to those observed in reflections), the attenuation is uniform as the number of layers is increased. 
These properties assure an almost diffraction-free propagation through the structure, with a similar 
attenuation and the size of PSF for symmetric and non-symmetric composition and for a broad 
range of structure sizes L. This behaviour is illustrated in fig. [3fci with an FDTD [16] simulation 
showing the uniform non-diverging distribution of the Poynting vector inside the structure for a 
beam size of the order of A/ 10. Transmission takes place along the direction normal to the layer 
boundaries with negligible divergence. Furthermore, our material may be used to fabricate slabs cut 
at arbitrary angle to layer surfaces. As is shown in fig. 03, c we continue to observe the diffraction- 
free propagation for inclined layers, whether the aperture covers only a single or multiple periods 
of the grating. This may be heuristically explained with the similar PSF for a symmetric and non- 
symmetric composition of multilayers. Both of them are encountered at the cross sections drawn 
along the propagating beam, starting from various points inside the aperture and normal to the 
layer boundaries. Furthermore, it is possible to construct other simple optical nano-elements. In 
fig. [H we demonstrate the operation of a double multilayer, which may be part of a cloaking device 
or an optical interconnect, as well as prism for imaging subwavelength beams. 
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4. CONCLUSIONS 

We have demonstrated a diffract ion- free, low-loss material, which is impedance matched to air. 
The diffraction-free propagation is verified by the analysis of FWHM of PSF in the function of the 
number of periods. The material consists of silver and a dielectric with refractive index of n — 2.6 
and operates at A = 422nm. It has the effective imaginary part of refractive index smaller than 
that of silver by a factor of 50. The reflections are in between — lOdB and — 30gLB, and the FWHM 
of PSF is at the order of 40nm. Small reflections, small attenuation, and reduced Fabry Perot 
resonances make it a flexible metamaterial for arbitrarily shaped optical planar devices with sizes 
of the order of one wavelength, such as the elements of optical interconnects or cloaks. 
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